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Derivative structural polytypes coexisting with the rhombohedral A7 structure of elemental bis¬ 
muth (Bi) have been discovered at ambient condition, based on microstructure analyses of pure Bi 
samples treated under high pressure and high temperature conditions. Three structures with atomic 
positions close to those of the A7 structure have been identified through first-principles calculations, 
showing these polytypes energetically comparable to the A7 structure under ambient condition. 

Simulated diffraction data are in excellent agreement with the experimental observations. We argue 
that previously reported variations in physical properties (e.g., density, melting point, electrical 
conductivity, and magnetism) in bismuth could be due to the formation of these polytypes. The 
coexistence of metastable derivative structural polytypes may be a widely occurring phenomenon 
in other elemental materials. 

PACS numbers: 81.30.Hd, 81.40.Vw, 61.50.Ah 

Bismuth as a typical semimetal lias attracted exten- earlier observations were attributed to a so-called “mo- 
sive research interests. Numerous studies have been con- saic block” structure [20, 21] with slightly different lattice 
ducted on Bi to investigate the electrical properties [1 constants [18]. A variety of formation mechanisms were 

7], thermoelectric properties [8, 9], and phase transitions proposed for this mosaic structure based on different as¬ 
under high pressure and high temperature (HP-HT) [10- sumptions, such as the influence of impurities, allotropic 

14], etc.. Under ambient condition, Bi is commonly des- changes, and deformation twinning [18, 22]. Nonethe- 

ignated as a rhombohedral lattice (space group R-3m, less, the atomic configuration of the mosaic structure in 

Strukturbericht A7), which is characterized by a pair of Bi crystals remains elusive. 

atoms spaced non-equidistantly along the trigonal axis j n this paper, we report structural polytypes of Bi 
in a Peierls distortion of the simple cubic structure [15]. evidenced by both experimental observations and first- 

Alternatively, the structure can be described either as principles calculations. We performed structural charac- 

a hexagonal lattice with six atoms per unit cell or as a terizations on pure Bi samples treated at various pres- 

pseudo-cubic lattice with one atom per unit cell [5]. How- sure and temperature conditions with XRD and high- 

ever, it is well-known that the observed X-ray diffraction resolution transmission electron microscopy (HRTEM) 

(XRD) patterns of Bi cannot be fully accounted for with an( j discovered a remarkable structural diversity in Bi 

the rhombohedral A7 structure, which is designated as after the treatments. Using crystal structure prediction 

(marked by) Doubtful (?) in the PDF card (JCPDS techniques, we identified three new structures of bismuth 

PDF ^44-1246) of the International Center for Diffrac- that are closely related to the well-known A7 structure 

tion Data and labeled with § in the periodic table of an d differ from each other only by minute lattice distor- 

elements [16]. The symbol of “§ indicates crystal struc- tions. First-principles calculations indicate these struc¬ 
ture is unusual or may require further explanation. tures have very similar energies, explaining the observed 

Some long-standing unusual behaviors in physical coexistence of multiple metastable polytypes in Bi. Such 

properties have also been noted in Bi, possibly result- a structural diversity is intrinsic, not due to extrinsic 

ing from the structure uncertainty. In 1928, Kapitza re- action from the impurities or deformation, and leads to 

ported inconsistent changes in specific resistance when Bi the property variations (such as density, melting point, 

crystals were subjected to strong magnetic field [17], and electrical conductivity, and magnetism) observed in Bi 

attributed the observations to small lattice distortions in ln ear li er studies. These polytypes provide a structural 

Bi structure. Later, slight density variations of Bi single interpretation for the so-called mosaic block structure 

crystals were reported [18]. In 1931, Goetz and Focke previously speculated. 

found that Bi single crystals grown under a strong ori- Bi granules (purity 99.999%, Alfa Aesar) were com- 
ented magnetic field exhibit a decrease in density and an pressed into cylinders (06mmx6mm ) and subjected to 
increase in specific resistance [19]. The applied magnetic HP-HT treatments in a six-anvil high pressure appa- 

field appeared to affect the Bi lattice during solidification, ratus [23]. The HP-HT conditions are presented in 

leading to changes in the density and resistance. All these [24](Table SI). After the HP-HT treatments, the samples 
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were quenched to room temperature before the pressure 
was unloaded. Both the raw granules and the bulk sam¬ 
ples after HP-HT treatment were ground into powders for 
XRD (Rigaku D/MAX-PC/2500, Cu-Ka) and transmis¬ 
sion electron microscopy (JEM-2010) characterizations. 

Searches for energetically and dynamically possible 
structures of Bi were carried out at pressures from 
0 to 6GPa with simulation cells consisting of 2-12 
atoms using CALYPSO, without employing any known 
structure information [25]. Structural relaxation was 
then performed using the density functional theory 
(DFT) within the generalized gradient approximation 
(GGA)[26]implemented in VASP [27]. Phonon frequen¬ 
cies were generated using the CASTEP with finite dis¬ 
placement theory [28]. The calculation parameters are 
provided in [24].It is known that conventional GGA cal¬ 
culations tend to overestimate the lattice parameters [29]; 
hence a hydrostatic pressure of 1.2 GPa was applied to 
the A7 structure during structural relation, producing 
lattice parameters consistent with the experimentally de¬ 
termined ones under ambient condition (PDF #44-1246). 
This pressure (1.2 GPa) was thus applied to all proposed 
structures in first-principles calculations. 

An HRTEM image of the sample treated at 2 GPa 
and 2273 K viewed along [010] zone axis is presented in 
Fig. la with the corresponding selected area electron 
diffraction (SAED) patterns shown in Fig. lb. Sim¬ 
ulated atomic distribution of the ideal A7-Bi structure 
viewed along same zone axis is shown in Fig. lc. In 
Fig. la, several regions with distinct atomic configura¬ 
tions are highlighted (yellow circles, labeled A-F), which 
are clearly different from the ideal A7 atomic distribu¬ 
tion (Fig. lc). The boundaries between these lattice 
distortion regions (Fig. la A F) are fuzzy and poorly 
defined. For clarity, these regions are displayed side-by- 
side in Fig. 2, and the d-spacings and interplanar angles 
measured from these regions are summarized in Table 
SII [24], All these regions show very similar d-spacings 
and interplanar angles, with deviations smaller than 0.1 
A and 0.3°, respectively. We note that the d-spacings ( 
3.74±0.03 A and 3.97±0.03 A ) and the interplanar an¬ 
gles ( 71.52±0.05° ) determined from region E perfectly 
match those of the (003) and (101) planes of the A7 struc¬ 
ture. Therefore, the SAED pattern (Fig. lb) taken from 
the entire area of Fig. la does not show clear evidence 
of structural variation, even though changes in lattice 
parameters and clearly differences in lattice fringes are 
observed. 

The shifts of Bi atoms from the A7 structure are so 
small ( less than 0.1 A ) that HRTEM image and the 
corresponding SAED patterns alone cannot offer solid 
proof of structural polytypes. To account for the sub¬ 
tle lattice distortions observed in HRTEM measurement, 
we theoretically investigated a number of possible struc¬ 
tural variations, in which atoms were shifted slightly from 
the equilibrium sites of the A7 structure without signifi¬ 



FIG. 1: TEM of Bi treated under 2GPa and 2273 K. (a) 
HRTEM image viewed along the [101] zone axis of the 
A7 structure. Different lattice distortion regions are 
labeled by A F. Measurements of d-spacings and 
interplanar angles for the A7 structure are shown in 
region E; (b) Overall SEAD pattern corresponding to 
(a); (c) Simulated atomic distribution observed along 
[010] direction in A7-Bi structure. 


FIG. 2: Enlarged atomic-resolution electron 
micrographs for A-F areas in HRTEM image shown in 
Fig.l. Corresponding atomic arrangements for given Bi 
structures are projected on the micrographs. 


cant energy increase. Three new structures based on CA¬ 
LYPSO searches along with the known A7 structure and 
phase II (structure of Bi at pressure higher than 2.7 GPa) 
are shown in Fig. 3, where each structure is depicted 
with biatomic layers to highlight the structural similari¬ 
ties among them. These newly identified structures (Fig. 
3b-d) are designated as M (monoclinic with space group 
C2/m), T1 (triclinic with space group P-1), and T2 (an¬ 
other triclinic structure with space group PI) phases, re¬ 
spectively. The symmetry groups, enthalpies, densities, 
and lattice parameters of these structures are listed in 
Table I. Unit cells and atomic positions are shown in 
Fig. SI and Table SIII[24], respectively. 
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FIG. 3: Top views of possible Bi polytype crystal 
structures, (a)-(e) are views along the a or b axis for 
A7, M, Tl, T2, and Bi-phase II, respectively. The 
double atomic layers are described by the aotms 
connevted by illustrative bonding sticks; (f) the typical 
chair-form liexatomic unit, and (g) the distorted 
chair-form hexatomic unit for the T2 structure. 


The bismuth biatomic layer can be constructed by 
chair-formed hex-atomic rings as building blocks. In Fig. 
3f, Bi atoms labelled as Bil, Bi3, and Bi5 in the hex- 
atomic ring constitute the upper sub-layer, while those 
labelled as Bi2, Bi4, and Bi6 constitute the lower sub¬ 
layer. The Bi structures are constructed by stacking a 
series of biatomic layers. The biatomic layers in A7, 
M, Tl, and Phase II structures are all parallel to the 
ab plane of the lattice. The A7 structure is comprised 
of • • ABCABC- • • stacked biatomic layers A, B, and C 
along the c axis, with each layer spanning one third of the 
lattice period. For M, Tl, and phase II structures, each 
biatomic layer spans one lattice period along the c axis. 
T2 structure (Fig. 3d), however, possesses the most com¬ 
plex crystal structure: The repeating unit consists of four 
chair-formed hex-atomic rings with distinct bond angles 
and bond lengths, as detailed in Fig. 3g. Each biatomic 
layer can only be constituted by moving the units along 
the [100] and [221] directions. These layers are parallel 
to the (024) plane of T2 structure, with an • • ABAB- • • 
stacking along the b axis. 

The chair-formed hex-atomic rings in these structures 
differ slightly in the bond angle and bond length. The 
corresponding angles (al-a6) and lengths (LI—L6) are 
thus introduced (see Fig. 3f), and listed in Table I to 


TABLE I: Calculated symmetry groups, lattice 
parameters, enthalpies, densities, and double-atomic 
layer characteristics of the structural polytypes. 


Structure 

Symmetry 

group 

lattice 

parameter 

Enthalpy 

(eV/atom 

Density 
) (g/cm3) 

layer structural 

characteristic 



a=b=4.5540 A 



o1=q2=q3=q4 

A7 

R-3M 

c=11.8281 A 

-3.8783 

9.82 

= o5 = o6=94.54° 


< d L> 

o=/3 = 90° 



L1 = L2 = L3=L4 


'7=90° 



= L5=L6=3.1 A 



a=7.8873 A 



ol = o2 = o4=o5 = 94.55° 

M 

C2/m 

b=4.5572 A 

-3.8777 

9.75 

o 3=o 6=94.70° 


< C 2D 

c=6.5836 A 



L1 = L4=3.102 A 


0 

o 

II 

i'¬ 

ll 

S 



L2=L3=L5 



13=143° 



=L6=3.098 A 



a=4.5738 A 



al=o4=94.98° 



b=4.5786 A 



o 2 = o 5=94.78° 

Tl 

P-1 

c=4.8088 A 

-3.8769 

9.47 

a3=o6=94.61° 


( o i) 

a=90.31 o 



L1 = L4=3.098 A 



/3 = 118.25° 



L2 = L5=3.123 A 



7=60.04° 



L3=L6=3.106 A 






o 1=94.92 ±0.07° 



a=4.5767 A 



o2 = 94.92±0.08° 



b=4.7786 A 



o4=94.91±0.06° 

T2 

PI 

c=13.3039 A 

-3.8781 

9.56 

o 5=94.88 ±0.07° 


(Cj) 

0=86.07° 



o3=o6=94.725 ±0.005° 



/3=90.02° 



L1 = L4=3.098 A 



7=90.02° 



L2=L5=3.1115±0.005 A 






L3=L6=3.110 A 



a=7.0015 A 



ol = o3=o4=o6 = 94.91° 

Phase II 

C2/m 

b=6.2556 A 

-3.8406 

10.07 

o2 = o5 = 145.28° 


< C 2 3 D 

c=3.3749 A 



L1 = L2 = L4 


0=7 = 90° 



=L5 = 3.277 A 



/3 = 112.71° 



L3=L6=3.199 A 


characterize each structure. The A7 structure has a 
unique chair conformation with identical bond angle of 
94.54° and bond length of 3.1 A . For the newly identified 
Bi polytypes, the bond angle and length show gradually 
increased diversities, leading to a gradually reduced sym¬ 
metry from M to Tl, and to T2. However, we note that 
deviations in the bond angle, bond length, and enthalpy 
are only within 0.5%, 0.8%, and 0.04%, respectively, for 
the Bi polytypes compared with the A7 structure. 

Dynamic stability of the new structures was verified 
by the phonon spectrum calculations (Fig. S2, [24]), 
which show no imaginary frequencies in the Brillouin 
zone. Thermodynamic stability of these derivative struc¬ 
tures was also examined, and the results are shown in 
Fig. 4. Enthalpy differences of the new structures rela¬ 
tive to the A7 structure are less than 3 meV at pressures 
between 0 and 2.7 GPa. Such small differences suggest 
all these new structures are stable within this pressure 
range. At pressures higher than 2.7 GPa, our calcula¬ 
tions show that the most stable structure is Bi-phase II. 
This is consistent with the known phase diagram of Bi 
[ 8 ], 

The very small enthalpy differences among the A7, M, 
Tl, and T2 structures suggest that very small energy 
perturbations, due to either phase transition from phase 
II on pressure release or magnetic treatments, may alter 
the structure of Bi and that these derivative polytypes 
can co-exist with the A7 structure. This may be the 
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FIG. 4: Enthalpy differences (AH) of the predicted Bi 
polytype structures relative to A7. Inset shows the 
details of AH at the pressure range from 2.2 to 3.2 GPa. 


structural origin of the speculated “mosaic block” or sec¬ 
ondary structure in pure Bi [18-21]. The density of calcu¬ 
lated A7 structure is 9.80 g/cm 3 , in agreement with the 
experimental value reported by Kapitza [17]. The den¬ 
sities of M, Tl, and T2 structures are all slightly lower 
than that of A7, with the lowest being 9.47 g/cm 3 in 
Tl structure. Thus, the coexistence of various structural 
polytypes explains well the density variations measured 
by Goetz in Bi single crystals solidified under different 
magnetic field strength [21]. 

Table SII summarizes d-spacings and interplanar an¬ 
gles measured from HRTEM for regions A-F of Fig. 2 
in comparison with the calculated ones [24]. Simulated 
atomic positions viewed along the corresponding direc¬ 
tions of the structural polytypes are superimposed on 
the HRTEM images in the upper right corners (Fig. 2). 
The excellent agreement between simulated and observed 
images suggest that region E corresponds to a projection 
along the [010] direction of the A7 structure, regions A 
and C to a projection along T2 [100], regions B and F to 
a projection along Tl [010], and region D to a projection 
along M [010]. 

The structural polytypes can also be detected in XRD 
patterns (Fig. 5), which further confirm the coexis¬ 
tence of the derivative polytypes in bulk samples. XRD 
peaks generally match the rhombohedral A7 structure 
(PDF$T4-1246), but with additional weak peaks near 
the main peaks that cannot be explained by the A7 struc¬ 
ture alone. The calculated XRD patterns of the M, Tl, 
T2, and A7 structures are also shown in Fig. 5 for com¬ 



20 30 40 50 60 70 80 26.8 27.2 


2-Theta (Deg.) 

FIG. 5: Observed XRD patterns of various Bi samples 
as compared with those calculated new from the new 
structural polytypes, (a) Experimental XRD spectra 
corresponding to the Bi samples obtained under raw 
state, cold-pressing at 2GPa and 2.7GPa, HT-HP at 
2GPa-1073 K and 2GPa-2273 K, the calculated XRD 
spectra of M, Tl, T2 and A7 structures, (b) The 
enlarged spectra at the range of 26.5-27.5°. 


parison. For illustration purposes, we concentrate on the 
main peak in the XRD patterns, which are magnified 
in Fig. 5b for the two-theta range from 26.5° to 27.5°. 
Clearly, there are one or more additional peaks other 
than the A7 (012) peak at 27.16° in all the samples. A 
small peak at 26.99° corresponding well to T2 (111), is 
present in the raw material (untreated) and the sample 
cold-press to 2 GPa, suggesting that small amounts of 
the derivative polytype T2 phase coexisted in the sam¬ 
ples. Two additinal peaks at 26.82° and 27.09° present in 
the the sample cold-press to 2.7 GPa (i.e., above the A7 
■fT phase II boundary of 2.52 GPa) can be attributed to 
the Tl (111) and M (112), respectively. In the samples 
treated at 2 GPa/1073 K and 2 GPa/2273 K, the peaks 
at 26.75° and 27.09° may correspond to the diffraction 
lines of Tl (100) and M (112), respectively (Fig. 5b). 
Note that the last two samples underwent melting and 
recrystallization, a more intense new peak is present at 
a smaller angle. 

Based on the experimental and theoretical investiga¬ 
tions, we establish the structural diversity as an essential 
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property for element Bi, in particular for samples treated 
at high pressures and high temperatures. Such diversity 
leads to the formation of secondary structures, which 
may be considered as derivative structures of the A7 
phase. Coexisting derivative structural polytypes may 
be a general phenomenon. One example is the mono¬ 
clinic diamond (M-diamond) recently reported in bulk 
nanotwinned diamond material [30]. In that case, phase 
transition from carbon onion nanoparticles to diamond 
and localized stress heterogeneity across diamond twin 
domains may be the cause. The exceptional thermal and 
mechanical properties observed in such multi-phase ma¬ 
terial suggest that coexisting structural polytypes may 
have interesting potentials in tailoring physical proper¬ 
ties of materials. 
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RESEARCH METHODS 


Experimental details 

Bi granules (Alfa Aesar, purity of 99.999%) were compressed into cylinders (6 mm in 
diameter and 6 mm in height) and each was placed into a BN capsule with 8 mm outer 
diameter. High-pressure and high-temperature (HP-HT) experiments were performed using 
a China-type large volume cubic press with a maximum of 1400 tons on every WC anvil [? ]. 
Pyrophyllite cubes with edge length of 49 mm served as both pressure medium and gasket. 
The samples were first compressed to required pressure, and then heated to the desired 
temperature for 1 hour by resistive heating with graphite tube heaters (outer diameter of 
10 mm and inner diameter of 8 mm). The HP-HT conditions are presented in Table SI. 
For all the experiments, temperature was measured in situ with a type C thermocouple 
(W5/Re26). Pressure was estimated with previously obtained calibration curves based on 
the electrical resistance change[? ] during the phase transitions of Bi at room temperature. 
After HP-HT treatment, the samples were quenched to room temperature at a cooling rate 
of ~100K/s before the pressure was unloaded. Bulk samples recovered at ambient condition 
had diameters of approximately 5.5 mm and heights of 5.3-5.7 mm. Both the raw granules 
and the bulk samples after HP-HT treatment were ground in an agate mortar with a pestle, 
for XRD (Rigaku D/MAX-PC/2500, operated at 40 kV and 200 rnA with a Cu-Ka target) 
and transmission electron microscopy (JEM-2010 with an accelerating voltage of 200 kV). 
The HR-TEM was carefully calibrated with silicon at the selected magnification. The change 
in microstructure introduced due to grinding is expected to be the same for all samples; thus 
the difference observed in TEM and XRD analyses are considered due to HPHT treatment. 

First-Principles Calculations 

The software package Crystal structure AnaLYsis by Particle Swarm Optimization (CA¬ 
LYPSO) [? ] was employed to search possible low-energy structures for Bi. We examined 
thousands of candidate structures with simulation cell sizes ranging from 2 to 12 atoms, at 
pressures of 0-6 GPa. Structural relaxations were performed using the density functional 
theory (DFT) within the generalized gradient approximation (GGA) implemented in the 
Vienna Ab Initio Simulation Package (VASP). Projector augmented-wave pseudopotential 
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with 6s 2 6p 3 electrons (as is the case of the valence for Bi) was adopted. The energy cutoff 
was 300 eV, and appropriate Monkhorst-Pack k meshes were chosen to ensure that en¬ 
thalpy calculations converged to better than 0.01 meV per atom for each structure. Phonon 
frequencies were generated using the Cambridge Serial Total Energy Package (CASTEP) 
with finite displacement theory. Conventional DFT schemes that apply to the GGA usually 
overestimate lattice constant parameters[? ]. In order to make direct comparison between 
calculated lattice parameters and the experimentally observed values, we adjusted hyro- 
static pressure in computing the lattice parameters of the A7 structure, until they match 
the parameters for the rhombohedral structure (PDF#44-1246). This hyrostatic pressure 
was identified for 1.2 GPa and then applied to all other possible structures in first-principles 
calculations. Finally, XRD patterns of these structures were simulated by Reflex tools pack¬ 
age of the Materials Studio (MS) software, to help interpreting the experimentally observed 
XRD patterns. 
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Table S I: Details of the HP-HT treatment conditions 


Sample No. 

Pressure (GPa) 

Temperature (K) 

Sample state under 
treatment conditions 

1 

24=0.1 

3034=5 

Phase I 

2 

2.74=0.1 

3034=5 

Phase II 

3 

24=0.1 

10734=10 

Liquid 

4 

24=0.1 

22734=30 

Liquid 


Table S II: Correspondence between the measured data for different areas in HR-TEM and 
corresponding d-spacings and angles in predicted Bi structures 


Measuring data in HR-TEM image 

Corresponding d-spacings and angles in 

predicted Bi structures 

Area in HR- 
TEM image 

d-spacing( A ) 
Deviation: 4=0.03 

Angle ( ° ) 
Deviation:=t0.05 

Structure and 

zone axis 

d-spacing( A ) 

Angle ( °) 

A 

di=3.99 

71.38 

T2-Bi 

d(oi 2)= 4 -0 

71.34 


d 2 = 3.76 


[100] 

d(oi2)=3-7453 


B 

d 1= 3.97 

71.54 

Tl-Bi 

d(ooi) = 4-0323 

71.50 


d 2 = 3.76 


[oio] 

d (10l) =3.7593 


C 

di=3.97 

71.36 

T2-Bi 

^(oi2)=4-0 

71.34 


d 2 = 3.76 


[100] 

^(oi2) = 4-7453 


D 

di=3.95 

71.67 

M-Bi 

d(ooi)=3.9683 

71.64 


d 2 = 3.76 


[010] 

^(20l)=3-7453 


E 

d 1= 3.97 

71.52 

A7-Bi 

d(003)=3-9370 

71.56 


d 2 = 3.74 


[010] 

d (10 i)=3.7453 


F 

di=3.97 

71.43 

Tl-Bi 

^(ooi)=4-0323 

71.50 


d 2 = 3.76 


[oio] 

d (10l) =3.7593 
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Table S III: Atomic positions in new structures of Bi 


Bi structures 

Atoms 

Wyckoff 

x/a 

Atomic Positions 

y/b 

z/c 

M 

Bil 

4i 

1.029(0) 

0 

0.793(1) 

T1 

Bil 

2i 

0.965(5) 

0.269(6) 

0.202(8) 


Bil 

la 

0.125(0) 

0.971(2) 

0.991(7) 


Bi2 

la 

0.624(6) 

0.528(2) 

0.508(5) 


Bi3 

la 

0.124(8) 

0.029(1) 

0.758(4) 

T2 

Bi4 

la 

0.625(1) 

0.470(7) 

0.741(8) 


Bi5 

la 

0.124(3) 

0.028(7) 

0.258(3) 


Bi6 

la 

0.624(6) 

0.470(4) 

0.241(9) 


Bi7 

la 

0.124(9) 

0.969(9) 

0.491(7) 


Bi8 

la 

0.624(6) 

0.529(5) 

0.008(5) 



Figure S 1: Unit cells of the Bi crystal structures. 
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-1 --- L 

G F Q Z 

T1 structure 



-1 ---- 

G F Q Z G 

T2 structure 



Figure S 2: Phonon spectra of new Bi structures. 
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